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CN Abstract 

<^^ In the last decade the subordinated processes have become popular and found many practical applications. Therefore 

04 in this paper we examine two processes related to time-changed (subordinated) classical Brownian motion with drift 

5—1 (called arithmetic Brownian motion). The first one, so called normal tempered stable, is related to the tempered 

^^ stable subordinator, while the second one - to the inverse tempered stable process. We compare the main properties 

^5 (such as probability density functions, Laplace transforms, ensemble averaged mean squared displacements) of such 

two subordinated processes and propose the parameters' estimation procedures. Moreover we calibrate the analyzed 

systems to real data related to indoor air quality. 
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1. Introduction 



Processes based on the Brownian motion were considered in many aspects and have found various practical ap- 
^-H plications [1-8]. But the assumption of normality for the observations seems not to be reasonable in the number of 

J> examined phenomenon. Therefore in many Gaussian models the Brownian motion is replaced by its various mod- 

^^ ifications. One of the simples modification is the extension of Gaussian by another distribution, for example the 

^^ Levy-stable one. Processes based on the stable distribution are very useful in modeling data that exhibit fat tails. 

(^— •) For example, the classical Omstein-Uhlenbeck process was extended to the stable case and analyzed in [9, 10] as a 

suitable model for financial data description, see also [11]. The another possibility of modification for Brownian-type 
processes is introduction a time-changed Brownian models. This extension is related to replacement of real time in 
^1 Brownian systems by non-decreasing Levy process (called subordinator), that in this case plays a role of random (op- 

T— I erational) time. The new process is called subordinate. The idea of subordination was introduced in 1949 by Bochner 

K^ [12] and expounded in his book in [13]. The theory of subordinated processes is also explored in details in [14]. 

The subordinated processes based on the diffusive Brownian motion were considered in many disciplines. The 
Brownian motion subordinated by gamma process, so called variance gamma, is analyzed in [15] in the context of 
3 option price modeling. The other applications and main characteristics of such system one can find in [16]. Subor- 

dination of Brownian motion by the inverse Gaussian process is called the normal inverse Gaussian (NIG) and was 
considered for example in [17] and proposed to modeling turbulence and financial data. Applications of NIG pro- 
cesses to asset returns are also shown in [18] while analysis of real environmental data by using NIG distribution is 
presented in [19]. Let us mention that the Brownian motion driven by strictly increasing Levy subordinator is also a 
Levy process, moreover when the subordinator is temporally homogeneous Markov then the subordinate has also this 
property, [14]. 

Another possibility of subordination is replacement the real time in Brownian difl'usion by inverse subordinators 
and processes that arise after this transformation are called anomalous diffusion. In the domain of anomalous diffusion 
the typical approach is based on continuous time random walk (CTRW), [20, 21], and subordinated Levy processes can 
be treat as a limit in distribution of CTRW, [22] . The key issue in the framework of CTRW as well as in subordination 
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technique is the waiting-times distribution corresponding to observed constant time periods [23]. In the last decade 
the anomalous diffusion processes were analyzed by various number of authors in many disciplines. For example 
the subordinated Brownian motion driven by inverse Levy-stable subordinator was considered in [11, 22, 24, 25], the 
inverse tempered stable subordinator was examined in [23, 26-28] while inverse gamma process was mentioned in 
[29]. The general case of Levy processes that can play a role of inverse subordinators were explored for example in 
[30-32]. 

In this paper we examine two processes related to subordinated Brownian motion. The first one, so called normal 
tempered stable [33, 34], is a Brownian motion with drift (called arithmetic Brownian motion, ABM) driven by 
tempered stable subordinator, while the second one is an ABM subordinated by inverse tempered stable subordinator. 
The tempered stable processes are extension of the a-stable Levy systems but possess also the properties of Gaussian 
models, therefore in the last few years they have become popular and very useful in description of many real data, 
[10, 34, 35]. We compare the main statistical properties of the ABM driven by tempered stable and inverse tempered 
stable subordinator. Moreover in two considered cases we propose the parameters' estimation procedures and validate 
them. In order to illustrate the theoretical results we calibrate the examined processes to real data related to indoor air 
quality. 

The rest of the paper is organized as follows: In section 2 we introduce the ABM and tempered stable subordina- 
tor. We present the main properties of those processes and define the time-changed ABM driven by tempered stable 
process. For this system we also examine the main statistical characteristics, such as probability density function, 
Laplace transform and ensemble averaged mean squared displacement that can be an useful tool to distinction be- 
tween diffusion and anomalous diffusion models. In this section we propose also the estimation procedure based on 
the distance between theoretical and empirical Laplace transforms and validate it. In section 3 we examine inverse 
tempered stable subordinator and its main statistical properties as well as define the ABM driven by inverse tempered 
stable process. In this section we also present the estimation procedure for unknown parameters. In section 4 we 
model real data sets related to indoor air quahty by using the mentioned subordinated processes. Last section contains 
conclusions. 

2. Arithmetic Brownian motion with tempered stable subordinator 

2.1. Arithmetic Brownian motion 

The arithmetic Brownian motion (ABM) is a process {X{t), f > 0) defined by [23, 36]: 

dX{t) = pdt + dB(t), (1) 

where {B(t), f > 0) is a classical Brownian motion and/3 > 0. The solution of equation (1) takes the form 

X(t) = X(0) + fit + B(t). 

In the further analysis we assume X{0) - with probability one and in this case the process defined above has 
Gaussian distribution with mean /3t and variance /. It is called also Brownian motion with drift, [23] and first of 
all was used to description of stock prices [1, 37]. It is an extension of the classical Brownian motion therefore its 
modifications have found many other practical applications, like diffusion in liquids modeling, [38] and description 
of hydrology time series [39]. 

2.2. Tempered stable subordinator 

The tempered stable subordinator {T{t), f > 0) is a strictly increasing Levy process with tempered stable incre- 
ments, i.e. with the following Laplace transform, [27, 28]: 

< g-^^W >= e'C-i'-C-i+^r)^ i > 0, < a < 1. (2) 

When A - Q, then {r(f)) becomes totally skewed a-stable Levy process. The probability density function (pdf) of 
tempered stable subordinator can be express in the following form: 
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where {U(t), f > 0) is a totally skewed a-stable Levy motion with the stability index a, [23, 40]. By using tail 
approximation of stable density, [41], we obtain the following: 

for some constant Ca- Therefore the right tail can be approximated by: 

where FT(t)(-) is the distribution function of T(t). Let us mention that for /I = (the a-stable case), the above formula 
reduces to 1 - Fij(,)(x) ~ x^" for fixed t. 

2.3. ABM with tempered stable subordinator 

The ABM driven by tempered stable subordinator is the process {Yrit), t >Q} defined as follows: 

Yrit) = X(T(t)), 

where {X{t)} is ABM defined in (1) and {r(f)) is strictly increasing Levy process with tempered stable increments 

with the Laplace transform given in (2). We assume those two processes are independent. 

Using the form of the solution of equation (1) we can write the exact formula for the process {Yx(t)}, namely: 

YT(t)^B(T(t))+/3T(t)- (3) 

The process {Yrit)} is known in the literature as normal tempered stable, [33] and the algorithm for generating its 
sample trajectories in points f i , t2, ...ytn proceeds follows: 

1 . Simulate the increments of subordinator dTi - T(ti) - T(ti-i) with the initial value Tq = 0. In order to do this 
use the simple algorithm presented for example in [27]. 

2. Simulate n standard Gaussian random variables A^i , N2, ■•■, N„ 

3. Define dYjii) = Ni y/dTi+fidTi and put Fr(f,) = ^Li ^^'■■ 

The sample trajectory of the process {Yrit)} with parameters a = 0.8, A = I and/3 = 1 is presented on the top panel of 
Fig.l. 

On the basis of the formula for Laplace transform of general subordinated processes given in [ 14] we can calculate 
this characteristic of the process {Yj(t)}, namely: 

< e-^i-HO >^ exp If I A"- -Ia+/3z- ^A ]| . (4) 

The pdf of the process [Yrit)} defined in (3) can be calculated by using the following formula: 



fYAt)(x) = fx(z)(x)fT(,)(z)dz, 

Jo 

where /x(z)(-) and fT(t)(') denote the pdfs of the ABM {X(t)} and tempered stable subordinator {r(f)), respectively. 
Therefore we obtain: 



poo 1 

fvMx)- ^e-^^-P^^'l^^-'^-'"'fu,,(z)dz, 

Jo ylnz 



ilnz 

where fu{t)i-) is the pdf of the totally skewed a-stable Levy motion {L'^(f), f > 0) mentioned above. 

On the basis of formulas (3) and (4) we can also calculate the main characteristics such as mean and autocovariance 

function: 

< Yrit) >^ Ptar-^ 
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Figure 1: The sample trajectory of the process {Yrit), t > 0] subordinated by tempered stable subordinator (top panel) and process {Ys{t), ^ > 0} 
subordinated by inverse tempered stable subordinator (bottom panel) with parameters a = 0.8, A = i 3nd/3 = 1 on the interval [0, 1]. 



cov(t, s) =< Yrit), Yt(s) >-< Yrit) >< Yt(s) >= min(s,t)[aA"-^ +/3^a(l - a)A"-^) . 

We consider also one of the most popular characteristic of the process that is especially important in real data analysis 
because it can be an useful tool for recognition between diffusion and anomalous diffusion models. This characteristic 
is called mean squared displacement (MSD). Here we consider the ensemble averaged MSD that is defined as a second 
moment of a given process. Therefore for {Yrit)} it takes the form, [29]: 



< Yj(t) >-- 



X fYj.(,)(x)dx, 



where fYj{t)(:) is a pdf of the process {YT(t)}. Let us emphasize, that the ensemble averaged MSD for real data can be 
calculated only in case when there are available more than one trajectory. For the process [Yjit)] defined in (3) the 
ensemble averaged MSD is a polynomial of the second order with respect to t: 



< Y^it) >= t\l3^aA"-^f + t[ar-^ +/3V(1 - a)A!'-^) . 



As we observe, in case of y6 = 0, the ensemble averaged MSD scales as t. In Fig. 2 (top panel) we present this 
characteristic calculated on the basis of 1000 trajectories for the ABM with tempered stable subordinator with a - 0.8, 
A - \ andjS = 0.01. Moreover we also show the fitted second-order polynomial. 
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Figure 2: The ensemble averaged MSD calculated on the basis of 1000 trajectories for the {Yrit)] (top panel) and {Ys (0) (bottom panel) processes 
withff = 0.8, A = 1 and/S = 0.01. On the top panel we present the fitted second-order polynomial while on the bottom panel - fitted power functions 
for small and large values. 



2.4. Estimation procedure 

We propose to estimate the parameters a, A and p by using the formula of the Laplace transform given in (4). 
Because [Yjit)] is a Levy process, [14], therefore its increments are independent and have the same distribution with 
the following characteristic: 



< e 



-z{YT(t+d)-YT(t)) 



>^&^^{d\r -{a+pz 



1 



The estimation scheme for random sample y\,y2, ■■■,yn proceeds as follows: 

1. Calculate increments of the observations: dyt = yi+i -y, for / = 1,2, ...,n - 1. 

2. For the increments find the empirical Laplace transform by using the following formula: 

^ «-i 

n - 1 ^—^ 

1=1 
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Figure 3: The boxplots of the values of estimators for parameters of ABM with tempered stable subordinator {Yrif)]. The values were calculated 
on the basis of 1000 trajectories of length 1000 each. The theoretical values are : a = 0.26, /i = 6, yS = 0. 11. 



3. By using the least squares method find the estimators a, A, p that satisfy: 



(a, A,I3) = mina,Aj3 U>(z) - exp U" - U + /3z - -z' 

In order to show the efliciency of the described method, in Fig. 3 we present the values of the estimated parameters 
for simulated process {Fj-(f)). To the analysis we take 1000 trajectories of length 1000. The theoretical values are: 
a - 0.26, A - 6, p - Q.W. As we observe the theoretical values are close to the estimated values (are between 
appropriate quantiles) that indicates the procedure works properly. 



3. ABM with inverse tempered stable subordinator 

In this section we examine the main characteristics of the ABM with driven by the inverse tempered stable subor- 
dinator. Such kind of systems were considered in the literature in many aspects, [23, 26-29], therefore some properties 
we present without proofs. 
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Figure 4: The sample trajectories of the tempered stable subordinator {T{t)] (top panel) and inverse tempered stable subordinator {S(t)| (bottom 
panel) for a = 0.26 and A = 6. 



3.1. Inverse tempered stable subordinator 

The inverse tempered stable subordinator is the process {5(t), t > 0) defined as follows: 

5(T) = inf{f >0: r(f) > t), 

where {T(t), f > 0) is the tempered stable Levy process defined via its Laplace transform in (2). In Fig.4 we present 
the trajectories of the tempered stable subordinator {r(f)) and corresponding tempered stable inverse subordinator 
{5(t)) for parameters a - 0.26 and /I = 6. As we observe, the processes are non-decreasing but exhibit completely 
different behavior. 

Using the relation P{S{t) < x) - P{T{x) > t) we can calculate the pdf of the process {5(t)), namely: 



fi 



d r 



Therefore we get: 



dx Jo 



where {U(t)} is a totally skewed a-stable Levy process described in the previous section. 
Because fu(x)iu) = -^fumiu/x^'"), [42], therefore we get: 

1 pT/x^'" pT/x'l" 

f^x Jo Jo 



c/m 



ax 

ffJc Jo " ^ Jo 



Jo 



As a final result we obtain: 



ax 
In case of /I = we have: 



/S(T)W = TfT(x)(T) + A ufT(x)(u)du - A" I fT(x)(u)du. (5) 

o'-*^ Jo Jo 



fs(T)(x) = Tfu(x)(T), 

ax 



that coincides with the result presented in [42] for a-stable case. For large t the density given in (5) tends to: 

fs(r)(x) ~ — T/r(,)(T) + ^— < T(x) > -A" = — T/rw(T). 
ax ax ax 

Similar, for small t we get: 

fs(T)(x) TfT(x)iT)- 

ax 
On the basis of equation (5) we can calculate also the Laplace transform of {5(t)), namely: 

Jo 
However such derivation requires numerical approximations. 

3.2. ABM driven by inverse tempered stable subordinator 

The ABM driven by the inverse tempered stable subordinator is defined as follows: 

Ys{t)^X(S{t))^I3S(t) + B{S(t)), 

where {X{t)] is the classical ABM explored in section 2. The simulation procedure of the process {Fs(t)), t > 0} 
is completely described in [27]. The sample trajectory of the process is presented in Fig. 1 (bottom panel). To the 
simulation we take the following values of the parameters: a - 0.8, A - \ and/? = L Let us remind that on the top 
panel we present the sample trajectory of the ABM with tempered stable subordinator {T{t)}. For {}s(t)) we observe 
constant time periods that is typical for subdiffusive processes. This behavior is not visible for the [Yjit)]. 
The Laplace transform of the ABM driven by inverse tempered stable subordinator [S (t)) can be calculated by using 
the following formula: 



Jo 



< e 



where fs(T)i') is given in (5). Thepdf /j'j(t)(x) of the process {Ys{t)} satisfies the generahzed Fokker-Planck equation: 

dfY,(T)(x) 



dr 



d 1 d^ 
dx 2 d^x 



^(/ys(T)(x)). 



where O is an operator defined in [23]. On the other hand we can use the formula for pdf of subordinated processes 
given in [14]: 

fY,(T)(x) = fx(z)(x)fs(T)(z)dz. 

Jo 

Therefore we obtain the following: 

r°° 1 

Jo ^27:2 
The mean of the process {Fs(t)) is given by, [29]: 

<Ys(t)>=/3<S(t)>=/3 f e--^"u''-^Ea,a((AuT)du, 

Jo 

where Eafiiz) is a generalized Mittag-Leffler function defined as follows, [43]: 

For small t, < Fs(t) >~ t" , while for large values of t the mean tends to r, [26]. Moreover the ensemble averaged 
MSD takes the form, [45]: 

< y|(T) >=/32 < s^{t) > + < B^(S(t)) >^/3^ < S^(t) > + <S(t)> 

where 

<S^(t)>= I X^fs(r)(x)dx. 

Jo 
For small t the last characteristic behaves: 

^Q CFX ^0 (XX Jo X 

On the other hand for large r we have: 

r*co 

Jo 
Therefore as a final result we obtain that the ensemble averaged MSD for the process {Ys{t)} satisfies: 

< Y^(t) >~ ciT^" + C2t'", when t -* 



< y^W >~ dir-'^e--"^ x-^'e" ^^i''> fumix)dx + diT, when t ^ oo 
Jo 

for some constants c\, C2, d\, d2- Let us mention ci and di depend on /3 in such way that for/3 = 0, ci = di = 0, 
therefore for /? = the ensemble averaged MSD behaves like t" for small values of t and like t for large values of 
this parameter, [26]. 

The ensemble averaged MSD calculated on the basis of 1000 trajectories for the process {1s(t), t > 0) we present 
on the bottom panel of Fig. 2. To the simulation we take a = 0.8, A - I and /3 - 0.01. We show also the fitted 
power function for small and large values of arguments. Let us mention that on the top panel we present the ensemble 
averaged MSD for the ABM with tempered stable subordinator {Yx(t), f > 0) with fitted polynomial of order 2. 
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Figure 5: The boxplots of the vakies of estimators for parameters of ABM with inverse tempered stable subordinator {Ys(t)]. The values were 
calculated on the basis of 1000 trajectories of length 1000 each. The theoretical values are : a = 0.4, A = 0.2, fi = 0. 



3.3. Estimation procedure 

The estimation procedure for parameters a, /3 and A of ABM with inverse tempered stable subordinator is partially 
described in [23]. We only mention here that it is based on the decomposition of the time series into two vectors. 
The first one is responsible for the lengths of constant time periods typical for the subdiffusive processes, while the 
second one appears after removing the constant periods. According to the theory, in our case the first vector constitutes 
independent identically distributed sample from tempered stable distribution. Because the right tail of the tempered 
stable distribution behaves like e^^'^x^" therefore we propose to estimate the parameters a and A by fitting this function 
to the empirical tail by using the least squares method. The second vector is related to the external process {X{t)}, 
therefore in the case of ABM the differenced series has Gaussian distribution with mean /?. Therefore the parameter 
we estimate as a mean of the differenced series. 

In order to present the efficiency of the presented method we simulate 1000 trajectories (each of length 1000) of ABM 
with inverse tempered stable subordinator {Ys{t)} with a = 0.4, A - 0.2 and/? = 0. Next, we estimate the parameters 
by using the presented scheme and in Fig. 5 we show the results. As we observe the fitted parameters are close to the 
theoretical values. 
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Figure 6: The humidity (top panel) and temperature (bottom panel) measured by three sensors between 9:48 and 23:59 on 12.08.2010. 



4. Applications 



In this section we examine real data sets that describe humidity (in %) and temperature (in °C) of the indoor air in 
some open space of huge company. For simplicity we denote humidity as DATAl and temperature - as DATA2. Those 
two analyzed quantities were measured by three sensors placed in the open space. Therefore for DATAl and DATA2 
we have three paths (trajectories). The humidity and temperature was quoted per minute. To the analysis we take data 
from 9:48-23:59 on 12.08.2010 (851 observations of each trajectory). On Fig. 6 we present paths of humidity (top 
panel) and data related to temperature (bottom panel). 

As we observe on Fig. 6, DATA2 exhibit behavior typical to subdiffusive process {^^(t)) described in previous 
section, namely the constant time periods. This behavior is not visible for DATAl (top panel). In the first step we 
examine the empirical ensemble averaged MSD calculated on the basis of three trajectories (of DATAl and DATA2), 
see Fig. 7. To the ensemble averaged MSD of DATAl we fit the second-order polynomial (top panel), while for small 
and large values of the MSD of DATA2 - the power functions. In both cases we use the least squares method. As 
we observe, the ensemble averaged MSD for DATAl behaves like second-order polynomial for all arguments that 
suggests behavior similar to this observed in Fig. 2 (top panel). The empirical MSD of the temperature that we 
observe on bottom panel of Fig. 7 clearly suggests subdiffusive behavior. 

In the next step of our analysis we estimate the parameters of the suggested processes, i.e. ABM with tempered 
stable subordinator for DATAl and ABM driven by inverse tempered stable subordinator for DATA2. In both cases 
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Figure 7: The ensemble averaged MSD of humidity (top panel) and temperature (bottom panel) together with the fitted second-order polynomial 
(for DATAl) and power function for small and large arguments (for DATA2). 



we use the procedures presented in sections 2 and 3, respectively. In the Table 1 we present the estimated parameters 
for DATAl. 



DATAl 


a 


A 


yS 


trajectory 1 


0.23 


6.1 


0.11 


trajectory 2 


0.29 


5.9 


0.12 


trajectory 3 


0.24 


6.1 


0.12 



Table 1: Estimated parameters of the ABM with tempered stable subordinator for DATAl. 



Moreover in Fig. 8 we present the empirical Laplace transform for each of three trajectories of DATAl and theoretical 
one given in (4). In order to confirm that the tempered stable subordinator is better than the stable one we show also 
the Laplace transform for the process {Yrit)} with A = (the other parameters we estimate by using the same method 
as this presented in section 2). Moreover we test also the Gaussian and a-stable distributions for differenced series 
of DATAl by using the methods based on the distance between empirical nad theoretical distribution functions, [44]. 
Those tests reject the hypothesis of the stable or Gaussian behavior of the observed data. 
After analysis of the humidity we start the examination of temperature. The ensemble averaged MSD indicates the 
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Figure 8: The empirical Laplace transform (solid line), theoretical Laplace transform for ABM with tempered stable subordinator (star line) and 
theoretical Laplace transform for ABM with stable subordinator (dashed line) for three trajectories of DATAl. 



data can be considered as subdiffusive process. Therefore we propose to use the ABM with inverse tempered stable 
subordinator described in section 3. Let us mention that the a-stable subordinator is not appropriate in this case 
because of the behavior of ensemble averaged MSD that in stable case behaves like t" for all values of arguments. 
The estimated values of a, A and /? parameters for three trajectories we present in Table 2. 



DATA2 


a 


X 


/S 


trajectory 1 


0.41 


0.13 


-0.0085 


trajectory 2 


0.3936 


0.11 


-0.0039 


trajectory 3 


0.4906 


0.129 


-0.012 



Table 2: Estimated parameters of the ABM with inverse tempered stable subordinator for DATA2. 



5. Conclusions 

In this paper we have examined two processes related to subordinated ABM driven by tempered stable type 
systems. The first one, so called normal tempered stable, arises after subordination of ABM by strictly increasing 
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tempered stable process, while the second one is a result of subordination of ABM with inverse tempered stable 
model. We have compared the main characteristics of such systems, like Laplace transforms, asymptotic behavior of 
pdf as well as the ensemble averaged MSD, that can be an useful tool of recognition between considered models. We 
have described also the estimation procedures for the parameters of considered processes and validate them. Finally, 
we have analyzed the real data related to indoor air quality in context of presented methodology. 
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